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ABSTRACT 

Investigations of interacting and merging galaxies at high redshift are vital to our 
understanding of their formation and evolution. To date, the identification of interac- 
tions at z ~ 3 and above has relied on rest-frame ultraviolet morphological parameters. 
Here, we present five serendipitous spectroscopic z ~ 3 Lyman break galaxy (LBG) 
pairs with projected proper separations < 15 /i~^kpc in our survey of nine separate 
deep Keck fields. The data consist of 140 of our highest signal-to-noise ratio LBG spec- 
tra and ~ 500 of our most confident colour-selected LBGs. We show that the pairs 
are composed of two distinct close and/or interacting LBGs from a detailed analysis 
of the rest-frame ultraviolet 1-D and 2-D spectra and the deep broadband images. 
In addition, we show that the number and separation distribution of the pairs is ex- 
pected from (1) the two-point angular correlation function when applied to the LBG 
pair separation distribution in our survey and ~ 2500 colour-selected LBGs from the 
literature and (2) an analysis of a carefully matched high-resolution hybrid numerical 
and analytical cosmological simulation. Because the spectroscopic slitlets have ran- 
dom orientations with respect to the close pairs on the sky, the serendipitous pairs 
provide an unbiased sampling of the underlying close pair fraction. Finally, we discover 
two Lya emitters (LAEs) in our slitlets and find that they reside within 50 projected 
/i~^kpc of the spectroscopic LBGs. In this work, we uncover a strong relationship be- 
tween Lya emission and pair separation. All confirmed and all candidate LBG pairs 
with separations of ^ 15 projected /i~^kpc exhibit Lya in emission and we find an 
indication of an overabundance of Lya emission in pairs with ^ 50 projected /i~^kpc 
separations. This relationship suggests a picture in which a measurable fraction of the 
Lya emission of LBGs, and potentially LAEs, is generated via interaction mechanisms 
such as triggered star formation and the dispersal of obscuring gas and dust. As a re- 
sult, serendipitous spectroscopic close pairs provide a unique means to help identify 
and study high-redshift galaxy interactions using ground-based optical data. 

Key words: methods: observational - galaxies: evolution - galaxies: formation - 
galaxies: high-redshift - galaxies: interactions 



1 INTRODUCTION 

A fundamental prediction of LCDM models of structure for- 
mation is the hierarchical build up of gala xies via m ergers 
with a rate that declines over time (e.g., Cole et aT]|200 0. and 
references therein). Observations of this process at low and 
high redshift are necessary to test these predictions, assess 
the effects of interactions, and to help clarify our current 
understanding of galaxy assembly. Close interactions and 
mergers act to trigger star formation and contribute to the 
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accretion of matter and the morphological transformation of 
galaxies over time (e.g., Larson & Tinsley 1978: 'Keel e t al.l 
1985; Mihos & Hernquist 1996; Barton et al. 2000, 20031 : 
ISteinmetz Navarrol 120021 : IScannapieco k Tisseral 120031 . 
and many others). 

In practice, however, identifying and studying merg- 
ers and interactions is difficult, largely because definitions 
vary, tidal features and interaction indicators can be faint, 
and morphologies are difficult to quantify. At higher red- 
shifts (z > 1), this is more challenging as galaxies become 
fainter and less resolved. Although nonparametric meth- 
ods have been developed to characterise the morphology 
of high-redshift systems (|Conselice et al.l l2003l : iLotz et aP 
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I2OO6I : iLaw et aP l2007l ) , such analyses are best done using 
the superior efficiency and resolution of optical data but, as 
a result, contribute an additional complication of linking the 
observed clumpy and chaotic rest-frame ultraviolet (UV) to 
the better understood rest-frame optical morphology. 

The Lyman break galaxies (LBGs) det ected by the Ly - 
man break colour-selection technique (Ste idel et alJ Il996l ) 
are a well- studied population of galaxies at high-redshift. 
Deep Hubble Space Telescope (HST) images of 2; > 2 LBGs 
in surveys such as GOODS, COSMOS, and the UDF show 
a wide range of rest-frame UV morphologies - from sin- 
gle compact or diffuse systems to multiple compact star 
forming components with, or witho ut, diffuse emission (e.g, 
iLaw et al.ll2007l : lRafelski et al.ll2009l ). The latter morphology 
has been interpreted as either star forming regions within a 
single coalescing galaxy and/or the result of interactions or 
mergers. 

Because of the faint nature of z 3 LBGs (L* corre- 
sponds to vsiR ^ 24.5), low- resolution spectroscopy can only 
achieve a signal-to- noise ratio (S/N) of a few with 8m-class 
telescopes for galaxies with m^ ^ 25.5 using reasonable 
integration times. This limitation prohibits detailed inves- 
tigations of individual LBGs. Except for the very brightest 
cases, detailed studies require composite spectra and/or rare 
gravit ationally lensed syst ems. Composite spectral analysis 
(e.g., "Shaplev et al.' '2003') has provided great insight into 
the average properties of LBGs but dilutes the features and 
relationships particular to individual systems. Gravitation- 
ally lensed LBG systems provide magnified flux that boosts 
the S/N of individual systems, but the total number of these 
systems is small, their reconstruction analyses are subject to 
model constraints, and their chance discoveries do not allow 
a study at desired redshifts. 

Recently, infrared integral field observations using adap- 
tive optics have been very informative on the kinematics of 
a sample of z ^ 2 LBGs (e.g., iForster Schreiber et al.ll200' 
iLaw et all l2007l : IShapiro et all I2OO8I : IWright et all lioq 
from prominent rest-frame optical features. Although this 
approach has been effective on relatively luminous LBGs, it 
would be prohibitively time consuming to compile a statisti- 
cal sample of typical (including faint) systems using existing 
instruments. Moreover, this method becomes increasingly 
difficult for higher redshift populations as key features are 
shifted out of the range of wavelength sensitivity. As a result 
of the above limitations, a clear and consistent interpreta- 
tion of the morphological and spectroscopic effect of mergers 
and interactions of typical LBGs at 2; > 3 via conventional 
approaches lay beyond the sensitivity thresholds and capa- 
bilities of current ground-based AO-assisted or space-based 
observations. 

A standard method to identify a large number of merger 
and interacting galaxy candidates at low- to intermediate 
redshift and to avoid many of the difficulties in morphologi- 
cal classification and int erpretation is to measure t he galaxy 
close pair fraction (e.g., [Patton et al.l 200C, I2OO2I: iLin et al ' 



20041 : iBerrier et al.1 12OO6I : iLin et al.1 l2007l : (Kartaltepe et al 
20071 ). Galaxy pairs with observed separations of ^ 50 
h~^kpc and velocity offsets of ±500 km s~^ are likely to in- 
teract and/or merge and can be more easily and consistently 
observed and quantified. Moreover, galaxies found in close 
pairs at high redshift have a higher likelihood of an imminent 
merger than at low redshift based in part on the fact that 



there are fewer galaxy clusters fe.g. jBerrier et aDl2006h . As 
a result, close pairs at high redshift provide a means to iden- 
tify galaxies that are interacting or are destined to interact 
on a short timescale and a strong constraint on the merger 
rate. Finally, the detailed behaviour of close pairs can be 
tracked in matched high resolution cosmological simulations 
to glean information beyond that of the observations. 

Although the expected number of LBGs in close pairs 
can be inferred photometrically from the observed angu- 
lar correlation function, spectroscopy is necessary to help 
differentiate actual close pairs from false pairs that have 
small projected separations on the sky. In order to gather 
a large number of LBG spectra, nearly all surveys to 
date have utilised optical low-resolution multi-object spec- 
troscopy (MOS). This method uses slitmasks to place slitlets 
at the locations of selected targets but cannot acquire spec- 
troscopy for objects with very small projected separations 
on the sky or that conflict in the dispersion direction. While 
MOS is efficient in acquiring a statistical sample of typi- 
cal LBGs, this shortcoming makes it difficult or impossible 
for a single slitmask to acquire information for many of the 
galaxies closely spaced on the sky. A programme employing 
a number of overlapping slitmask observations of the same 
area of the sky could, in principle, resolve this issue. How- 
ever, we show below that not all close galaxy pairs at z ^ 3 
are easily identifiable from their imaging alone and would 
introduce systematics in such a time expensive endeavour. 

A different technique, which we introduce here, allows 
existing MOS datasets to be used to study close and inter- 
acting pairs at high redshift. The density and clustering be- 
haviour of LBGs predict that a small fraction of close pairs 
are expected to fall serendipitously into the area probed by 
the MOS slitlets. Moreover, because the slitlets in conven- 
tional LBG surveys have random orientations with respect 
to the serendipitous close pairs on the sky, they produce 
samples unbiased by selection. In this manner, the serendip- 
itous pairs reflect the underlying close and interacting LBG 
pair behaviour and effectively circumvent the MOS mechan- 
ical constraints. 

We analyse a conventionally acquired nine-field deep 
Keck imaging and spectroscopic survey ioi z ^ 3 LBGs and 
find five serendipitous spectroscopic close pairs. We use 509 
colour-selected and 140 spectroscopically conf irmed LBGs, 
the large z ^ 3 photometric LBG dataset of Stei del et al.l 
(2003), and a high-resolution cosmological simulation to in- 
vestigate the distribution of close pair separations and the 
expectation of the serendipitous pairs. We present several 
lines of evidence from a detailed analysis of the imaging and 
spectroscopic data that support the close pair/interacting 
nature of the serendipitous detections. We show that the 
properties of this unique sample, in turn, provide insight 
into the rest-frame UV properties of the LBG population 
that are otherwise inaccessible and a means to help identify 
high-redshift interactions in ground-based data. 

This paper is organised as follows: We describe rele- 
vant information regarding the data and data acquisition in 
321 We present the serendipitous spectroscopic interacting 
pairs, investigate their rest-frame UV behaviour, and com- 
pare the results to precedence in the literature in ^Sl Close 
pair observations and measurements are discussed in 311 and 
the results from the simulation analysis in 33 The two Lya 
emitters discovered in our sample are presented in ^Sl and 
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we address the observed relationship between Lya in emis- 
sion and close pair separation in 33 Finally in ^ we provide 
a su mmary. All magnitudes are in the AB (Fukugita et al. 
1 199 6) magnitude system unless otherwise noted. We assume 
an Qm = 0.3, Ha = 0.7, /i = 1.0 cosmology throughout and 
list close pair separations in physical (proper) units. 



2 OBSERVATIONS 

We use the z ^ 3 LBG survey of ICooke et al. 1 (|2005l . here- 
after C05) for the observational analysis in this paper. 
The survey consists of deep imaging and spectroscopy of 
nine widely-separated fields covering ~ 460 arcmin^. Seven 
fields were imag ed using the Low Resolution Imaging Sp ec- 
trometer fLRIS: lOke et al.lll995l : iMcCarthv et al. 1998) on 
the 10-m Keck I telescope with two fields imaged using 
the Carnegie Observatories S pectrograph and Multi- Object 
Imaging Camera fCOSMIC: iKells et al.lll998l l mounted on 
the 5-m Hale telescope at the Palomar Observatory. We ob- 
tained follow-up spectroscopy of all nine fields using LRIS. 
The multiple fields suppress the effects of cosmic variance 
and provide ^ 800 colour-selected and ~ 200 spectroscopi- 
cally identified LBGs for statistical analysis. Details of the 
observations , data reduction, an d LBG analysis can be found 
in C05 and Coo ke et al.l (|2006l ) . Below, we briefly discuss a 
few relevant specifics of the data. 

The z ^ 3 LBG candidates in our survey are selected 
by their i^^BVRI colours. Although the deep images can se- 
lect LBGs to fainter magnitudes, we restrict our analysis 
in this paper to LBGs with m^ ^ 25.5 in order to define 
a sample that has (1) photometry accurate to ^ 0.2 mag 
uncertainty in all filters that probe the rest-frame UV con- 
tinua, (2) follow-up spectroscopy, (3) defined photometric 
selection functions, and (4) a measured spatial correlation 
function. 

The spectroscopic component of our survey follows a 
conventional approach to gather a large number of z ^ 3 
LBG spectra using MOS slitmasks. The survey was not 
originally designed to study interacting LBGs or to mea- 
sure close spectroscopic pairs. The physical constraints of 
MOS slitmasks can prevent acquisition of galaxies that are 
closely spaced on the sky or those that fall in the disper- 
sion direction of the targeted observations. This situation is 
illustrated in Figure [1] for two spectroscopically confirmed 
LBGs and associated colour-selected close pair candidates 
from our survey. However, because each MOS slit let covers 
a non- negligible area of the sky 18 arcsec^ on average), 
occasionally a close pair will have a small enough separation 
or coincidental near alignment to be detected in a slitlet that 
was originally intended for a single photometrically identi- 
fied colour- selected LBG candidate. 

The MOS slitlets in our survey were centred on the LBG 
candidates with position angles near the parallactic angle of 
the observations to minimise wavelength-dependent flux loss 
caused by atmospheric dispersion. As a result, this is an un- 
biased spectroscopic survey for the purposes discussed here 
because the data were acquired randomly with respect to 
the orientation of close pairs on the sky. Any detected close 
pairs in the slitlets are indeed serendipitous and directly re- 
flect the true underlying distribution of close pairs on the 
scales probed by the dimensions of our slitlets. 




Figure 1. Images of spectroscopic LBG close pair candidates 
(~ 20^' on a side) centred on a spectroscopically confirmed z ~ 3 
galaxy Left: LBG 1057-3275, labelled 3275, at z = 3.014 with 
colour-selected close pair galaxies 1057-3278 and 1057-3327 cir- 
cled and labelled similarly. Right: LBG 1643-961 at ^ = 3.412 
with the colour-selected close pair candidates 1643-979, 1643-949, 
and 1643-922. This figure illustrates how the physical constraints 
of multi-object slitmasks can make it difficult to acquire spec- 
troscopy of close galaxy pairs. The central rectangle replicates 
the spectroscopic slitlet dimensions and orientation for that par- 
ticular observation. The larger outer circle represents 50 h~^kpc 
physical radius from the central galaxy. The region bounded by 
the lines that extend from the upper and lower edges of the slitlet 
indicate the dispersion direction that is inaccessible for that ob- 
servation. For these particular observations, we were unable to ob- 
tain spectroscopy for candidates 1643-949, 1057-3278, and 1057- 
3327. Other constraints, such as the allowance of space between 
slitlets necessary for milling purposes, sufficient slitlet lengths for 
proper sky subtraction either side of the targeted object, and the 
presence of higher priority candidates in the dispersion direction 
(determined in the original survey) prevented the spectroscopic 
acquisition of close pair candidates such as 1643-922. 



The seeing in the images used for colour selection of the 
LBG sample ranged from 0.^^6 — 1.^^5 arcsec FWHM, with 
those having 0.''6 - 1.''0 FWHM used for pair analysis. In 
addition, all spectroscopy was acquired under good condi- 
tions (0.''8 - 1.''2 arcsec seeing FWHM). The slitlets were 
milled to have l.^^O — 1.^^5 widths, with minimum and 
^ lA" average lengths. With knowledge that z ^ 3 LBGs 
are near-point sources in ground-based observations, we de- 
signed the widths of our MOS slitlets 0.^^4 larger than the 
expected seeing FWHM to gather the maximum flux and to 
compensate for ^."2 — {)."3 astrometric errors across the 
focal plane. As a result, the spectral resolution for a given 
LBG was typically determined by the seeing and is ^ 400 
km s~^ FWHM for most observations. 



3 SERENDIPITOUS SPECTROSCOPIC CLOSE 
PAIRS 

In the full spectroscopic dataset, we discover LBGs that ex- 
hibit unexpected double Lya emission in the 2-D spectra 
with double or overlapping continua. These data form dis- 
tinct double Lya emission peaks in the merged 1-D spectra. 
In addition, we find that these objects show double morphol- 
ogy (two measurably separated flux peaks) in the broadband 
images. Furthermore, we find that the Lya emission and 
continua in the 2-D spectra have a one-to-one correspon- 
dence in the spatial direction when compared to the separa- 
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tions of the flux peaks in the images. We term these objects 
"serendipitous spectroscopic close pairs" (SSCpairs). 

A systematic analysis the images, 2-D spectra, and 1- 
D spectra of 140 of our highest signal- to- noise ratio (S/N) 
spectroscopically confirmed z ^ 3 LBGs finds five such sys- 
tems that show compelling evidence for more than one iden- 
tified major component. We discuss the observed proper- 
ties of the five SSCpairs in the following sections. Because 
the separations for some of the SSCpairs cause them to 
be unresolved by the so urce extraction routine SExtractor 
(jBertin k Arnoutd Il996') in the ground-based images, we 
adopt a naming convention that uses the R.A. of the tar- 
geted field and the identification number assigned by SEx- 
tractor for the originally targeted LBG. For the SSCpairs 
that were detected as two distinct LBGs by SExtractor, we 
label the pair using the R.A. of the field and the identi- 
fication number of the targeted LBG only. Figure [2] dis- 
plays the ground-based LRIS R-band images and 2-D spec- 
tra of these systems. Object 1643-2377 received follow-up 
higher S/N, higher reso lution longslit spectroscopy using 
LRIS (ICooke et al.l l20Q8V These data provide valuable in- 
sight into LBG close pair properties and are used below to 
help confirm our analysis methods of the low- resolution, low 
S/N SSCpair spectra. For the results presented here, includ- 
ing those from the simulation and close pair analysis, we use 
the parameters of the low-resolution data of 1643-2377 for 
consistency. 

In addition to the five SSCpairs, we find six LBG close 
pair candidates with separations < 15 /i~^kpc. These were 
not included as SSCpairs by definition because the bulk of 
the emission of the secondary galaxy did not fall into the 
dimensions of the slitlet. Four of the six SSCpair candidates 
are shown in Figure (3] with the remaining two shown in 
Figure [1] Future spectroscopy will determine whether the 
six candidates are indeed associated. 

The small separations and overlap of the faint SSCpair 
spectra make the clean extraction of individual coroponents 
difficult. As a result, we present the merged spectrgy in Fig- 
ure m along with spectral segments focusing on the Lya 
features. Presented in this way, the two components of each 
system can be directly compared, but have blended noise, 
continua, and features that include absorption effects from 
the Lya forest blueward of Lya emission. Each component 
of the SSCpairs shows Lya in emission and we show in ^3.21 
that this is not a selection bias. The presence of Lya emission 
from each component results in double-peak Lya emission 
profiles for all five merged SSCpair spectra. 

We remark that the double-peak Lya emission features 
result from merging the spectra and differ from the double- 
peak Lya feature predicted for a stati c or expanding shell 
of ga s arising from a single LBG (e.g., Tenor io-Tagle et al.l 



Il999l : iMas-Hesse et al.1 [20031 : IVerhamme et al 



2006 



2008 h. 

We describe below that the Lya emission features of the 
SSCpairs do not result from single galactic-scale expand- 
ing shells as is evidenced by the two spatially offset Lya 

^ Pair 0957-0941 has the largest separation of the SSCpairs. We 
made efforts to extract the individual components of the low S /N 
spectra using conventional routines, but found each spectrum was 
contaminated with a non-negligible amount of flux from the other 
galaxy. We present the merged spectra for this pair in Figure [H 
for consistency. 
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Figure 2. Ground-based R-band images and 2-D spectra of the 
flve spectroscopic z ~ 3 LBG close pairs, the SSCpairs. Each 
image is ~ 15^^ on a side and the labelled rectangle replicates 
the slitlet dimensions and orientation for that particular obser- 
vation. Directly below each image is a section of the correspond- 
ing 2-D spectrum centred on the Lyo; feature. Top: SSCpairs 
0056-0635 (0056-0635 and 0056-0632) and 0336-1250. Centre: 
SSCpairs 1013-1182 and 1643-2377. Bottom: SSCpairs 0957-941 
(0957-0941 and 0957-0961). The SSCpairs are named after the 
original targeted LBG. Objects 0056-0632 and 0957-0961 are de- 
tected as separate colour-selected LBGs but are referred to as 
0056-0635 and 0957-0941 for consistency. R-band images probe 
z ~ 3 LBG UV continua near rest-frame ~ 1700A. Although 
faint, the 2-D spectrum of 0056-0635 shows distinct continua and 
Lya emission for objects 0056-0635 and 0056-0632. 
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Figure 3. Ground-based R-band images and 2-D spectra of four 
z ~ 3 LBG < 15 h~^kpc close pair candidates. These four can- 
didates, and the two candidates shown in Figure ^ were not 
included in the SSCpairs analysis because the bulk of the flux 
from the secondary component fell outside of the slitlet dimen- 
sions. Top: Objects 0336-0380 and 1643-1687. Bottom: Objects 
2342-2525 and 1057-2462. The photometric close pair candidates 
are circled and labelled in the images. Object 1057-2469 is two 
objects, the fainter of which is the photometric candidate. As in 
Figure [2] each image is ~ 15^^ on a side and the labelled rectangle 
replicates the slitlet dimensions and orientation for that partic- 
ular observation. The images were obtained using LRIS except 
object 2342-2525 which was obtained using COSMIC. The four 
LBGs here, and two in Figure^ display dominant Lyo; emission. 



emission peaks and continua in each SSCpair 2-D spectrum, 
double morphology, and the direct correspondence between 
the spatial offsets of the 2-D Lya features and continua and 
the flux peaks in the images. We note that an accurate mea- 
surement of the small Lya emission velocity offset in the 
1-D spectra of SSCpair 1643-2377 is difficult. In that case, 
we used the high-S/N, higher resolution spectroscopic data 
(Figure [5]) to cleanly measure the Lya peak offset and find 
that the measurement agrees with our low-resolution value 
within the uncertainty. 

Table [1] lists the relevant information of the five SS- 
Cpairs. All reported velocity differences between the two 
components are corrected for offsets in wavelength resulting 
from their spatial separation as determined by the associ- 
ated fiux peaks in the images (§ I3.2|) . We see that the two 
components have projected physical separations of '-^ 3 — 13 
h~^kpc in the images and 2-D spectra. We remark that be- 
cause these are projected separations, they represent the 
minimum actual separations. 

O and B stars are largely responsible for the rest-frame 
UV flux over the wavelengths probed by the optical broad- 
band filters used in this survey. As a result, we interpret 
the continua and emission lines from the each SSCpair com- 
ponent as a major young star forming region or a distinct 
LBG. In the following sections, we discuss whether the re- 



sults favour an interpretation that the SSCpairs are coalesc- 
ing single massive galaxies with two violent starbursts or 
an interpretation that the SSCpairs are two separate close 
and/or interacting galaxies. The lines of evidence presented 
below argue that these are separate components that are 
inconsistent with the single galaxy interpretation. We use 
morphological and spectroscopic information obtained from 
the SSCpairs, space-based investigations of local LBG ana- 
logues, and LBG ground- and space-based observations and 
relationships found in the literature to help interpret the 
SSCpairs and to help shed light on the rest-frame UV be- 
haviour of LBGs. 



3.1 Morphology 

Space-based images of z - 3 LBGs (e.g., iLaw et al.|[2007h 
show a wide range of morphologies, such as concentrated 
single nucleated sources, multiple nucleated sources, and dif- 
fuse structures. Althou gh complex, typical LBGs have half- 
liRht radii of < 0.''3 ( Gardner et al.ll2000l : iFerguson etHI 
I2OO4I I , or < 2 /i~^kpc, and appear as near-point sources in 
ground-based images. We inspect the fiux contours of the 
SSCpairs in the LRIS ground-based i^^BVRI imaged for de- 
partures from an unresolved source in the form of significant 
extended features, elongations, or objects with more than 
one discernible fiux peak. 

The VRI broadband filters probe the star forming rest- 
frame continuum of the z ^ 3 LBG sample longward of Lya. 
We find that the centroids of the fiux peaks in the ground- 
based VRI images are accurate to ~ 0.^^4 (approximately the 
seeing half-width half maximum) aided by the near point- 
source profiles of LBGs. Each LBG that meets the 2-D and 
1-D criteria discussed below exhibits two measurably sepa- 
rated fiux peaks in the VRI images and have projected cen- 
troid separations of ^ 3 — 13 /i~^kpc. The separations for 
SSCpairs 1013-1182 and 1643-2377 are at or near the reso- 
lution limit. The double nature of 1013-1182 relies largely 
on other measurements presented in this work such as the 
two spatially offset Lya emission peaks in the 2-D spectra. 
For 1643-2377, the high S/N from the intrinsic luminosity of 
this system allows a secure centroid measurement. Finally, 
we note that the broadband continuum luminosity and Lya 
emission from each SSCpair component falls in the range 
of typical m^ ^ 25.5 LBGs at z ^ 3. Consequently, each 
SSCpair component has LBG-equivalent estimated star for- 
mation rates 10^ - lO^M© yr"\ Papovich et al. 200ll; 
IShaDlev et al.ll200lL l2003h . 

At low to intermediate redshifts, the UV and op- 
t ica l studies of compact UV- luminous galaxies (UVLGs; 
e.g., Heckman et al. II2QO5I : iGrimes et al]l2QQ7l : lOverzier et aD 
2008, hereafter O08) have made great strides in our under- 
standing the observed LBG UV morphology at high red- 
shift. UVLGs are in every practical way local analogues to 



^ LBGs at z ~ 3 are colour-selected as w^-band "drop-outs" be- 
cause the decrement in their continua shortward of the Lyman 
limit (912A, rest-frame) is redshifted into the w^-band. Therefore 
z ~ 3 LBGs are extremely faint (R ^ 27) or not detected in this 
bandpass. Inspection of the n^-band images reinforces the non- 
detection of both components detected in the longer-wavelength 
bands and helps to eliminate low-redshift line-of-sight interlopers. 
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Table 1. Serendipitous spectroscopic close pairs 



Pair 


R.A." 


DEC." 




^Lya 


ZISM^ 






Slit^ 




Seeing-^ 


0056-0635 


00:59:27.722 


+01:41:37.37 


25.0, 25.2 


3.367,3.386 


3.363,3.374 


870 


350 


1.5 


9.5 


0.8; 


1.4 


0336-1250 


03:39:17.199 


-01:31:29.10 


23.4^ 


2.824,2.845 


2.819,2.835 


960 


575 


1.5 


8.2 


0.9-1. 


1; 0.6 


0957-0941 


09:57:48.414 


+33:11:28.43 


24.8, 25.0 


2.680,2.688 


2.678,2.685 


480 


400 


1.5 


13.3 


0.9; 


1.0 


1013-1182 


10:15:57.234 


+00:21:55.06 


24.6^ 


2.779,2.792 


2.775,2.786 


730 


555 


1.3 


2.6 


0.9; 


0.9 


1643-2377^ 


16:44:48.314 


+46:27:08.23 


22.2^ 


3.036,3.044 


3.029,3.036 


430 


360 


1.5 


4.7 


0.8; 


0.7 



" Coordinates of the LBG flux centroid used to miU the multi-object spectroscopic slitlet. 

^ Assigned redshifts for the respective features of the two identifled components 

^ Velocity offset in km s"-*^ between the respective features, corrected for angular separation 

^ Width of multi-object spectroscopic slitlet in arcsec 

^ Projected angular component separation in the images in physical kpc 

f Spectroscopic FWHM in arcsec derived from night sky emission lines; R-band seeing FWHM in arcsec 

^ Integrated magnitude of both sources 

^ Values are from the low-resolution data for consistency 



LBGs having similar mass, gas content, half-light radii, UV 
colours, metallicity, and star formation rates. O08 find that 
although the objects in their sample do not show convincing 
evidence for mergers or interactions from the HST rest-frame 
UV imaging alone, every UVLG in their sample shows evi- 
dence of interaction in the rest-frame optical images. They 
conclude that mergers are the main mechanism for the ob- 
served star formation in these systems and find that the bulk 
of the star formation arises in several compact 100 — 300 
pc) regions. 

O08 deconvolve the pixel scale and dim the surface 
brightness of their z ^ 0.1 UVLG sample to match the 
depth and resolution of z ^ 3 LBG observations in the 
GOODS, COSMOS, and UDF fields. The multiple bright 
star forming regions of the UVLGs blend into 1—3 bright 
clumps with separations of < 0.5^' when redshifted to z 3 
and accurately reproduce the space-based morphology of 
LBGs. For the ground-based image resolution in our survey, 
these clumps would appear as a single near-point source for 
all but the most separated cases, which would be detected 
as SSCpairs depending on component properties. Perhaps 
the most convincing evidence of this is the case of SDSS 
J080844.26+394852.4 (hereafter UVLG 0808) in the O08 
sample. 

We find that all UVLGs in the O08 sample, except for 
UVLG 0808, have compact rest- frame UV morphologies and, 
if placed at z ^ 3, would be detected as single near-point 
source^ in our survey. UVLG 0808 has a fainter compan- 
ion ~ 2" to the SE that exhibits star formation flux in the 
UV image. An equivalent separation at z ^ 3 would be dis- 
cernible in ground-based images. In addition, the companion 
has a compact strong detection in the rest-frame Ha image 
and is the only UVLG in the sample that is detected as two 
separate galaxies by the Sloan Digital Sky Survey (SDSS). 

We search the SDSS spectra of the O08 sample for ev- 

^ Object SDSS J092600.41+442736.1 shows rest-frame Ha emis- 
sion with a separation that is at the limit of what is discernible 
as two objects in our ground-based survey. It did not receive HST 
rest-frame UV imaging, so it is unknown at this point whether it 
has two-component UV flux that would be detectable as a close 
pair at z ~ 3. Strong detection in the rest-frame UV is unlikely 
because one component is a much redder galaxy that shows an 
old, or dusty, population in the SDSS spectrum. 



idence of double-peak emission in the optical nebular lines. 
These lines are not subject to outflows as are Lya emission 
Unes f ^3.2.1[) and should closely trace the systemic veloci- 
ties. The multiple star forming regions of seven of the eight 
systems are located within ^ 1'' and the component rela- 
tive velocities are near or below the spectral resolution of 
SDSS, thus these UVLGs show single-peak emission. How- 
ever, UVLG 0808 shows evidence of a double-peak feature as 
witnessed in the Balmer lines. Because the two components 
of UVLG 0808 have an ^ 2" separation, a portion of rest- 
frame optical flux from the companion galaxy is expected to 
fall into the SDSS fibre and result in the weaker offset emis- 
sion line peaks {Av ^ 280 km s~^). UVLG 0808 is the only 
object in the sample that (1) has a well-separated compan- 
ion, (2) appears as two star forming point source galaxies in 
the rest-frame UV, and (3) shows evidence of double-peak 
emission in the 1-D spectroscopy. Observed at z ^ 3, a sys- 
tem with similar properties would appear by our definition 
as an SSCpair. 

Projected component separations of > 4 h~^kpc are 
necessary for individual identification of UVLGs in the O08 
sample and LBGs in our survey at the resolution of SDSS 
and LRIS, respectively. The remarkable similarities between 
the properties of UVLGs at low to intermediate redshift and 
LBGs at high redshift, and the evidence that interactions 
provide the bulk of the star formation in UVLGs, leads us 
to conclude that the separations and emission profiles ob- 
served in our SSCpairs likely result from two interacting 
LBGs. We investigate this further in the following sections 
and discuss the spectroscopic and morphological properties 
of LBGs found in the literature. 

3.2 Spectroscopy 

The 2-D spectra are individually examined for signs of mul- 
tiple spectra beyond the usual pipeline reduction of the tar- 
geted objects and distinct serendipitous objects. We search 
for evidence of closely spaced or overlapping spectra, and/or 
multiple Lya emission or absorption features. Alongside the 
image and 2-D analysis, we inspect the 1-D spectra for ev- 
idence of more then one system from redshift fits to ^ 20 
expected interstellar features. The actual number of features 
fit for a given LBG depends on the redshift and wavelength 
range probed by the MOS slitlet and does not include lines 
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Figure 4. Left: Spectra of the five LBG serendipitous spectroscopic close pairs shifted to the rest-frame and box-car smoothed by 3 
pixels. As discussed in the text, these spectra are the merged 1-D spectra of the two very close or overlapping spectra. Although this aids 
in the identification and comparison of the systems, it blends and can dilute the significance of individual galaxy features. Dashed (blue) 
and dotted (maroon) vertical lines mark the prominent ISM transitions of the apparent lower and higher redshift galaxy, respectively, 
and are labelled. For clarity, not all lines used in the fit are shown. The red (grey) line refiects the Icr per pixel error array, the green 
(light grey) horizontal line indicates zero fiux, and the thick solid (grey) vertical lines mark the positions of bright night sky emission 
lines. Right: 1-D spectral segments of the five pairs (not smoothed) that focus on a 24A region (rest-frame) centred around the merged 
double-peak Lya emission features. 



that fall near bright night sky emission features which can 
be difficult to extract cleanly from the faint spectra. The 
spectra in our sample have continua with relatively low S/N 
('^ 3 — 10) and are typical of LBG observations at z ^ 3. 

Lya is the most prominent rest-frame UV feature in 
z ^ 3 LBGs and, in contrast, can be detected at a high sig- 
nificance. The LBG spectra compiled to date show roughly 
50% of all LBGs dominated by Lya in emis sion, with the 
remaining dominated by Lya in absorption ("Shaple v et al.l 
I2OO 3, and this sample). Net EW values range from strong 
emission, > lOOA EW, to damped absorption (column den- 
sities of N(Hi) ^ 10^°-^ atoms cm"^). Here we term LBGs 
that are dominated by Lya emission in their spectra, eLBGs, 
and those dominated by Lya in absorption, aLBGs. After a 
thorough search for SSCpairs of all types, we observe that 
each SSCpair component, with no exception, displays Lya 
in emission. Thus, in our nomenclature, all SSCpair compo- 
nents are eLBGs. 

The Lya emission features and continua in the 2-D spec- 



tra for every SSCpair have the same separations in the spa- 
tial direction as the image centroid separation along that 
same calculated direction. The one-to-one correspondence 
helps to confirm the distinct two-component nature of the 
SSCpairs and enables Lya emission and star-forming con- 
tinua assignment to each peak in the images. We use the im- 
age centroids when correcting the Lya and ISM absorption- 
line velocity offsets for their angular separation. The separa- 
tions in the dispersion direction are complicated by the fact 
that Lya emission has a large range of velocity offsets with 
respect to the systemic redshift (see ^ I3.2.ip . Nevertheless, in 
the cases where the flux peaks in the images have small sep- 
arations in the spatial direction, large differences between 
the offsets of the two Lya features can provide additional 
evidence for two major star-forming components. 

In practice, detecting double Lya emission can prove to 
be easier than double Lya absorption when component ve- 
locity offsets and angular separations are small. Thus, it is 
possible, that even after a thorough search, we are missing a 
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Figure 5. Spectra comparing the Lyo; emission for 1643-2377 
obtained with two different resolutions. Left: The merged low- 
resolution discovery spectrum (~ 9 A resolution). Right: The 
higher signal-to-noise higher-resolution (~ 4 A resolution) follow- 
up spectrum showing two distinct peaks. 



fraction of LBG spectroscopic close pairs that only exhibit 
Lya in absorption because of the low S/N of the spectra 
and broad nature of the Lya absorption features. However, 
LBGs displaying Lya in absorption typically have stronger 
ISM transitions that can facilitate the de tection of a velocity 
offset between two overlapping spectra (jShaplev et al.ll2003l . 
C05). We find no aLBGs at larger radii where both the an- 
gular separation and velocity offsets would resolve the Lya 
features of the two systems and no clear evidence for double 
sets of ISM absorption transitions redward of Lya in aLBG 
spectra. In addition, we do not find SSCpairs having one 
component in emission and one in absorption in either the 
1-D or 2-D spectra. Finally, all spectroscopic components of 
the SSCpair candidates (Figure [T]&[3| show dominant Lya 
emission, with none exhibiting dominant Lya absorption. 
These findings are discussed further in JTl As a result, the 
observation that all SSCpairs are eLBGs is not primarily a 
selection effect. 



3.2.1 Lya offsets 

We measure corrected Lya velocity offsets for the SSC- 
pairs of AvLya 400 — 1000 km s~^. To expand on the 
difficulties in rest-frame UV continuum interpretation men- 
tioned earlier, Lya emission in LBGs is observed to be red- 
shifted by ^ 450 ± 300 km s~^ with respect to the systemi c 
redshift (C05 |Adelberger et al. 2003; Shaul ev et al]l2003l l. 
The velocity offsets likely result from galactic-scal e outflows 
caused by stellar and superno va- driven winds fe.g.. lHeckmanl 
l2002l : lAdelberger et al.l l2003h. Lya photons are absorbed by 
the approaching outflow while Lya photons scattered off- 
resonance from the receding shell can escape. This geom- 
etry results in Lya emission peak velocity offsets that are 
redshifted with respect to the systemic redshifts and de- 
pendent on the host galaxy properties. Consequently, the 
observed difference between the peaks of the Lya features 



cannot be used directly as a measure of the systemic velocity 
offsets of the SSCpair components. 

However, under the following assumption, we can in- 
vestigate the nature of the component Lya offsets using the 
Lya offset distribution of the LBG population. Because of 
the one-to-one correspondence between the ^3 — 13 /i~^kpc 
projected separations of the flux peaks in the images and the 
Lya peaks in the 2-D spectra, we attribute the double-peak 
Lya features to two distinct galactic-scale outflows. Admit- 
tedly, this is a simplifled picture but the relatively coarse 
resolution of our ground-based observations blends the ve- 
locities of any localised regions and can only probe net global 
outflows. 

"Adelberger et alj (|2003l l flnd a velocity offset distribu- 
tion between Lya emission and the ISM absorption lines 
(AvLya-iSM) of AvLya-iSM = 614 di 316 km s~"^. Bccausc 
LBG ISM absorption lines are blueshifted on average with 
respect to the systemic redshift by ^ —150 =b 160 km s~^ 
fAdelberger et all l2003l : IShaplev et aD 12003', and this sur- 
vey), this results in the above-mentioned 450 =b 300 km s~^ 
Lya v elocity offset distribution. In addition, Ishaplev et al.l 
(|2003l ) find that LBGs with dominant Lya emission (Group 
3 and Group 4 in that sample, or eLBGs in our nomen- 
clature) show smaller values of AvLya-iSM and therefore 
evidence of weaker outfiows. If we assume that eLBGs have 
a similar velocity offset distribution in behaviour and form 
to that of the full LBG population, the eLBG-only distribu- 
tion becomes AvLya-iSM — 518 ± 266 km s"\ with eLBG 
Lya offsets from the systemic redshift of 406 it 262 km s~^. 

If the Lya emission in the SSCpairs refiects outfiow from 
two star forming regions in single coalescing systems, we 
would expect the velocity difference between the two Lya 
features, AvLya, to randomly sample either of the above 
Lya velocity offset distributions. A random sampling of two 
values using the distribution of the eLBGs and full LBG pop- 
ulation results in expected differences of AvLya = 295 ±220 
km s~^ and AvLycx =~ 340 ±255 km s~^, respectively. How- 
ever, we find a much larger AvLya = 694 ± 234 km s~"^ 
for the SSCpairs. If we add a line-of-sight velocity compo- 
nent originating from two regions in a bound, massive LBG 
system by pulling from a < 250 k m s~^ distribution of rela- 
tive velocities (e .g. , Erb et al. 200i l2006l : ICenzel et aT]|2008l : 
iLaw et aI1l2009[ ) , we find an expected net two-component ve- 
locity offset distribution of AvLya ^ 320 ± 230 km s~^ (or 
~ 360 ± 260 km s~^ for the full LBG population), which 
falls short of describing our data. 

The correspondence between the 2-D Lya emission and 
continua and the two fiux peaks in the images forces us to 
conclude that the Lya emission features of the SSCpairs 
arise from either (1) single galaxies, but do not randomly 
sample the behaviour observed in the full or subset LBG 
population and that larger than typical AvLya values are 
correlated with star forming regions having substantial line- 
of-sight velocity contributions or (2) star forming regions 
with velocity offsets too large to be single systems and are 
truly pairs of LBGs. In either scenario, the double-peak Lya 
emission of the SSCpairs lends insight into the behaviour of 
z ^ 3 systems. 
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3.2.2 ISM absorption 

Spectra showing two distinct sets of ISM absorption hnes 
having large velocity offsets together with separate 2-D spa- 
tial offsets provide strong evidence for two independent 
systems. We cross-correlate and visually inspect the SS- 
Cpairs for ^20 strong ISM UV transitions that include 
Ly-a, Ly-/3, and strong ISM features such as SiiiA1260, 
Oi/SiiiAA1302,1304, CiiA1334, SiivAA1393, 1402, SiiiA1526, 
and ClvAA1548, 1550 A. We find evidence for two systemic 
redshifts from both the low- and high-ion transitions in the 
merged 1-D spectra. The ISM velocity offsets for the SS- 
Cpairs are listed in Table [T] We caution that the low S /N 
3—10) of the continua and the potentially complex nature 
of the rest-frame UV behaviour of the ISM, make these the 
least confident measurements in terms of significance when 
compared to the other analyses presented here. However, the 
ISM measurements of the higher S/N low-resolution spectra 
of 0336-1250 (- 6cr) and 1643-2377 (- lla) are more secure. 

Our data provide tests of the validity of these low S/N 
measurements. Object 0336-1250 was observed with three 
separate slit masks for a total integration time of 11700s and 
the luminous 1643-2377 was observed with three exposures 
of 1500s each. In an effort to mimic the lower S/N of the 
other SSCpair total observations, we performed fits to each 
of the three 0336-1250 slitmask observations and the individ- 
ual exposures of 1643-2377 separately. We found agreement 
to Az ^ 0.0015 between the fit to the individual exposures 
to their respective combined spectrum. 

The SSCpairs show a corrected velocity difference be- 
tween the two identified sets of ISM lines of Avism = 
448 lb 109 km s~^. As mentioned above, ISM lines exhibit 
an ^ 150 di 160 km s~^ range in blueshift from systemic. If 
we assume that the two sets of ISM lines randomly sample 
this distribution, we would expect Avism = 180 =b 135 km 
s~^. If we assume the SSCpairs were single systems with the 
bulk of the ISM absorption seen in the sightlines to two well- 
separated star forming regions with bound relative velocities 
< 250 km s~^, we would expect the combined effect to result 
in AvisM 194 ± 145 km s~^, which is much smaller than 
the observed value. As a result, the evidence in the ISM lines 
argue for the SSCpairs to be separate systems. 



3.2.3 Ly-a - ISM velocity difference 

A comparison of the difference between the SSCpair Lya 
and ISM features with the relationships found for LBGs 
is a further test of the viability of the 1-D spectral fea- 
tures to help confirm the dou ble-galaxy nat u re of the SS- 
Cpairs. As discussed above, I Shaplev et al.l (|2003l ) deter- 
mined AvLya-iSM = 614 zb 316 km s~^ for all LBGs and 
518 ± 266 km s"^ for eLBGs that are sim- 
ilar to the SSCpairs. Because we see two Lya emission 
peaks and evidence for two sets of ISM absorption lines, 
we treat the SSCpair components as individual systems 
and measure the AvLya-iSM for each component separately 
using no angular separation correction. This results in a 
AvLya-iSM — 475 db 205 km s""*^ when assigning the high- 
est redshift Lya feature to the highest redshift set of ISM 
features and the lowest redshift Lya feature to the lowest 
redshift ISM features. A reassignment of the Lya and ISM 
features results in the same central value but a larger scat- 



ter. The AvLya-iSM distribution found for the SSCpairs 
agrees with that of the LBG popul ation. Moreover, the dis- 
tribution follows the trend found in "Shaple v et aD (|2003h for 
LBGs displaying dominant Lya emission. 

In treating the SSCpairs as close pairs throughout 
this section, we would expect the 1-D spectral features of 
each galaxy to reflect the relationships and trends of the 
LBG population. We find that the AvLya, Avism, and 
AvLya-iSM are in good agreement with this picture and 
are at odds with an interpretation that the SSCpairs repre- 
sent starforming regions in single systems. 

3.3 Comparison to previous related work 

We look for observational precedent in the literature to help 
inter pret the S SCpairs by first exploring the LBG sample 
of S haplev et al. (2006). This sample consists of deep spec- 
troscopy of 14 LBGs at 2; ^ 3 gathered for the purpose 
of measuring the fraction of escaping continuum flux short- 
ward of the Lyman limit. The observations are valuable for 
an assessment here because they were selected in the same 
manner, lay within the same redshift path, and have spec- 
troscopy obtained using the same instrument (LRIS). More- 
over, the observations are beneficial here because they have 
a slightly higher resolution over the bulk of the more relevant 
wavelengths (< 1400A rest-frame) and have ^ lOx longer 
exposure times for S/N of ~ 5 — 25. In addition, the asso- 
ciated imaging data were obtained with similar or, for the 
subset of space-based data, superior depth and resolution. 

Of the 14 LBGs in that sample, nine are eLBGs. The 
Lya EW distribution of the sample is very similar to that of 
our full spectroscopic LBG sample (K-S test, p — 0.6). One 
object, D3, is considered a single system by I Shaplev et al.l 
(|2006 ) but consists of two galaxies (and separate spectra) 
with an angular separation of 1.^^9 (~ 15 /i~^kpc ). Three 
of the nine eLBGs are reported to show double-peak Lya 
emission. Because both components of D3 show Lya emis- 
sion, merging the two spectra would create a double-peak 
spectrum, making it a fourth double-peak system (similar 
to 0957-0941 in our sampl e). Relevant to later discussion, 
of the two LBGs shown by I Shaplev et al.l (|2006l l to exhibit 
measurable Lyman continuum flux, one is the double-peak 
Lya emission system C49 and the other is the pair D3. 
Narrow-band imaging of D3 shows that it is detected as 
a Lya "blob" in the sample of Matsuda et al.l (I2OO4) with 
Lya emission extending over 17 arcsec^. 

We now analyse the double-peak systems of 
IShaplev et al. I 12^3) as we did for the SSCpairs using 
the published information. Looking at the velocity offsets 
of the Lya feature of the four double-peak systems, we find 
AvLyc. = 533 =b 170 km s"\ cf. AvLyc. = 694 ± 234 km s"^ 
for our sample. Shapley and coworkers measure one ISM 
absorption redshift for each of the three Lya double-peak 
systems (pair D3 has separate values). If we assign the 
reported ISM absorption redshifts to the highest redshift 
Lya peak, we find AvLya-iSM = 1165 ± 667 km s~"^. This 
places these systems at the high-end tail of the AvLya-iSM 
distribution for the full LBG population, with one system 
having the highest (1823 km s~^) velocity difference of 
the entire z ^ 3 LBG spectroscopic sample. Assigning the 
features in this manner implies that it is highly improbable 
that the double-peak systems are single LBGs. 
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Instead, if we assign the ISM absorption redshifts to the 
lowest redshift Lya peaks we find AvLya-iSM = 579 =b 682 
km s~^. While the central value matches the expectations 
of eLBGs and is in good agreement with LBGs in gen- 
eral, the large scatter of the three systems is unrepresen- 
tative of the ±275 km s~^ scatter that we find for three- 
galaxy random samples taken from the LBG and eLBG 
data. In addition, we see similar relative offset and scat- 
ter behaviour when we assign the higher redshift ISM fea- 
tures to the lower redshift Lya emission (and vice-versa) 
in our SSCpair sample (section [3 .2.3 p . When inspecting the 
double-peak system C32, we find that the lowest redshift 
Lya peak is blueshifted with respect to the ISM absorp- 
tion lines by —183 km s~^. If we assign only this system 
to the higher Lya peak, then the three systems exhibit 
AvLycx-iSM = 803 =b 322 km s~"^ and the scatter is re- 
duced to the level of the expectations. When including the 
identified pair D3 {AvLya-iSM = 579 km s~^) in the Lya 
double-peak sample, we find AvLya-iSM = 713 =b 266 km 
s~^, which is more consistent with the distribution and scat- 
ter of the eLBG and LBG population. Although matching 
ISM features to one of the Lya emission peaks can make 
the double-peak systems fall in line with typical LBGs, what 
should we make of the second Lya feature? 

The above analysis shows that it is more likely that 
the three Lya double-peak LBGs are close pairs similar to 
D3 and the SSCpair sample as opposed to single systems. 
Such an interpretation produces a AvLya-iSM distribution 
in better agreement with the expectations of the LBG popu- 
lation and significantly reduces the scatter. It is possible that 
the presence of a se cond set of ISM absorption lines in the 
IShaplev et aD (|2006l ) sample may have been unexpected and 
ignored, identified and interpreted as strong ISM mixing, or 
identified and considered of too low significance as a result 
of blending. Although the S/N of the Shaplev et al. (2003) 
LBG sample is higher than what is typical for z ^ 3 LBG 
spectra, it does remain relatively low in that ISM absorption 
lines remain difficult to study with great confidence. Future 
investigation of the data (for example, by probing nebular 
lines in the IR) could resolve this matter. 

We remark that IShaplev et al.1 (|2003h find 33% (44% 
depending on slit orientation for D3) of the eLBGs in their 
sample, and 21% (29%) of the entire sample, to be double- 
peak systems that meet our SSCpair criteria. The data imply 
that higher S/N and/or higher resolution should reveal a 
larger fraction of SSCpairs than that seen serendipitously in 
our low-S/N, low-resolution survey, with the true fraction 
potentially even higher. 

Next, we examine the study of iLaw et al.l (|2007h that 
consists of space-based imaging of 216 z 2 — 3 LBGs and 
ground-based spectroscopy using LRIS with a similar con- 
figuration as our survey. In that work, the morphology of 
LBGs is quantified using a non-parametric analysis within 
1.5^^ radius (~ 13 kpc) which directly reflects the regime of 
the SSCpairs. In addition, this is the regime where systems 
are most sensitive to interaction that can visibly change their 
morphology and star formation properties. Law and cowork- 
ers investigate the rest-frame UV morphology of their LBG 
sample by quantifying the Gini parameter, colour dispersion, 
multiplicity, and size. The Gini parameter is a measure of 
the nebulosity or nucleated nature of the system, the colour 
dispersion is measured using V —z broadband colours and 



corresponds to rest-frame ^ 1500 — 2100A at z ^ 3, and 
the multiplicity parameter quantifies the multiple compo- 
nent riature_of the sources. 

iLaw et al.' (2007) investigate the relationships between 
the morphological parameters and the spectroscopic fea- 
tures Lya EW, low- and high-ionisation ISM EW, and 
ALy-a-iSM- As uotcd lu that work, only Lya EW shows 
a significant and consistent trend at z ^ 3 with all four 
morphological parameters. The Lya relationships indicate 
greater Lya emission EW with decreasing system size, de- 
creasing nebulosity (increasing Gini parameter), decreasing 
number of nuclei (multiplicity) or nuclei separation, and 
bluer colours. In addition, and in reference to low-redshift 
UVLGs, a visual inspection of the space-based images shows 
that systems with bright nuclei typically have one or two 
components whereas systems with faint nuclei typically have 
three or more components and are found in more diffuse sys- 
tems. 

Lya EW is seen to be strongest for large values of the 
Gini parameter indicating that nucleated sources are more 
likely to be eLBGs and diffuse sources are more likely to be 
aLBGs. The Lya EW increases monotonically over the mul- 
tiplicity range with values that indicate two strong double 
nucleated sources (multiplicity parameter value '-^10) down 
to a single nucleated source (multiplicity parameter value 
0). Finally, the Lya EW increases as the size of the systems 
decrease and the colour becomes bluer. iLaw et aP (|2007h 
show that each of these trends at z 3 has a significance 
of ^ 98%. Taken together, these data can be interpreted 
to show that double nucleated LBGs exhibit Lya emission 
and that the Lya emission EW increases as the components 
become bluer and as the separation between the two com- 
ponents decreases, or in other words, as the overall object 
size becomes progressively smaller for a given multiplicity. 

The SSCpair observations are consistent with this be- 
haviour and the interpretation that the SSCpairs are inter- 
acting bright double nuclei, or two single nuclei LBGs. The 
double nucleated sources with larger separations are those 
discernible as SSCpairs in ground-based data, as also in- 
dicated by the work of O08. Diffuse systems, and diffuse 
systems with multiple fainter clumps, would be observed as 
single LBGs in ground-based data. Moreover, such systems 
are LBGs with dominant Lya absorption from the Law et 
al. relationships. Strong Lya EW and blue continua can be 
evidence for recent star formation and/or decrease in dust 
and gas obscuration. The Lya EW of the SSCpairs in or- 
der of decreasing projected separation is 26, 31, 102, 56, and 
91 A. Although the SSCpairs sample is too small to deter- 
mine a trend, the SSCpairs have simil ar separations as the 
double nuclei LBGs in the Law et al ] (Hqoi) sample, have 
blue colours, and exhibit Lya emission with EWs that are 
in the upper quartile of the full LBG sa mple. 

We comment that the work of Law et al.l (|2007l l does 
not investigate the behaviour of LBGs with separations 
larger than ~ 13 kpc. We probe this regime in S] and ad- 
dress statistics applicable to this discussion in 33 In the next 
section, we test the expectations of finding the SSCpairs by 
measuring the distribution of LBG close pair separations in 
our full survey and a larger survey from the literature and by 
testing LCDM predictions using a high-resolution cosmolog- 
ical simulation. We show that the number of SSCpairs are 
indeed expected and, as a result, the investigation of this 
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and other similar samples offers a unique window into the 
behaviour and properties of close/interacting high redshift 
galaxies. 



4 SERENDIPITOUS PAIR EXPECTATIONS 
FROM CLOSE PAIR SEPARATIONS 

To investigate the distribution of LBGs on small scales (< 
500 /i-^kpc), we measure the number of LBG pairs versus 
separation from observation. To help test the universality 
of the results and to improve the statistics, we augment an 
investigation of our survey with the ph otometric and spatia l 
information from the larger survey of ISteidel et al. 1 (|2003), 
hereafter 803, consisting of ~ 2300 colour-selected and ~ 
800 spectroscopically identified LBGs in 17 separate fields. 



Cooke et al. (2005) dataset 




Separation (h 'kpc) 



4.1 Pair definition 

For the C05 survey, we consider only LBGs that are identi- 
fied in our images by SExtractor with uir ^ 25.5. The LBGs 
were selected using n^BVRI colour criteria 

{u-B)ab> 1-1 (1) 

0.6 < {B-R)ab < 2.1 (2) 

0.6 < {B-I)ab < 2.1 (3) 

{u - V)ab > 1.6 (4) 
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r Steidel et al. (2003) dataset 
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(V - R)ab < 0.6 (5) 

{V - I)ab < 0.6. (6) 

In order to provide more reliable photometric statistics, 
we use only colour-selected LBGs that meet the full crite- 
ria in the seven of the nine fields that have high-quality 
five filter imaging and LBGs that meet four-filter criteria in 
the remaining two fields (see C05 for details regarding the 
i^^BVRI selection criteria and the specifics of the color crite- 
ria for each of the nine fields). For the S03 survey, we use all 
LBGs that meet the {G-IZ) ^ 1.2, {Un -G)^ (G-IZ) + 1.0 
used in that work. These color-selection criteria are tested 
using > 300 and > 1300 spectra in the respective samples. 

LBG colour selection at z 3 is very efficient in filtering 
out objects not in the redshift path of interest (interlopers). 
As determined from the two spectroscopic samples, the frac- 
tion of interlopers in each survey is similar 0.20) when 
considering the particular treatment of the unidentified low 
S/N spectra. Differences in the photometric selection func- 
tions and selection efficiencies between the two surveys are 
small and are quantified and corrected for using the informa- 
tion from the spectroscopic samples. Figures [6] & [71 present 
the close pairs versus separation for this survey, the survey 
of S03, and the combined dataset. The two figures empha- 
sise large and small scales, respectively. The fall-off in the 
number of pairs at radii < 10 h~^kpc (best seen in Figure [T]) 
results from the resolution of the images and the ability for 
SExtractor to discern close galaxies as separate. 
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Figure 6. Histograms of pair separations versus radius. Top: 
Number of LBG pairs in our survey binned by 5 h~^kpc. Over- 
laid, and more distinguishable on small scales in Figure [Jl is a 
linear fit (thin solid red line), a power-law fit (thick dotted black 
line), and the corrected two-point angular correlation function of 
the form u;{6) = A6~^-^^ (thick solid blue line) fit to the data (see 
text). The contribution from LBG pairs to the angular correlation 
function fit is shown by the thick dashed (light blue) line, with the 
random interloper pairs comprising the remaining contribution. 
The thin (red) dashed line indicates the form of the random pair 
contribution normalised to data. As a comparison, the LBG con- 
tribution for an angular correlation function with slope 6 = 0.8 
(thick dot-dashed green line) is shown. Centre a nd bottom: 
The s ame as the top panel, but for the survey of ISteidel et al.l 



12 Cooke et al. 



4.2 Pair separations 

Overlaid on Figures [6] & [7] are the pair separation expec- 
tations from a linear (solid thin red line) and power-law 
(dotted black Une) fit to the raw data. We find that a power 
law slope of ^ = 0.53 ±0.01 provides a best fit for all three 
samples. As illustrated in the two figures, both forms fit well 
over all separations, with the linear fit showing a departure 
from the data near, and below, ^ 30 h~^kpc. Continuations 
of both fits argue for the existence of LBG pairs with sepa- 
rations smaller than the imaging resolution. 

As mentioned above, interlopers comprise a non- 
negligible fraction of colour-selected LBGs. To estimate the 
form of the contribution to the pair separations from random 
interloper line-of-sight projections and to estimate the like- 
lihood that the SSCpairs are random projections, we con- 
struct mock catalogues having the physical dimensions of 
the images. We insert galaxies with random locations in the 
plane of the sky following the number densities found for 
each field and the redshift distributions for each survey. 

Using catalogues of this design is effective at the small 
separations studied here because the likelihood of random 
pairs is low as a result of the relative low surface density of 
colour- selected galaxies (^ 1.7 arcmin"^). In addition, the 
enhancement of the number of random projected close pairs 
as a result of the spatial clustering of LBGs is negligible as 
compared to random Une-of-sight projections. The z ^ 3 
colour-selection techniques used here probe > 500 /i~^Mpc. 
The spatial correlation length of z ^ 3 LBGs is ^ 4 /i~^Mpc 
(|Adelberger et al.ll2003l . [20051 : ICooke et al.ll2005l ) with an av- 
erage of ^ 1.5 LBGs over random integrated over that vol- 
ume. Beyond that distance, the clustering effects of LBGs 
diminishes and the effects of random projections become 
dominant. Because the projected separations of the SSC- 
pairs are < 15 /i"^kpc , the increase in likelihood that the 
SSCpairs are random projections enhanced by the spatial 
clustering of LBGs is vanishingly small. 

We measure the interloper contribution for each field 
separately and combine the results from 1000 runs of the 
mock catalogues. The form of the contamination is shown on 
Figures[6]and[7] (dashed thin red line) normalised to the total 
number of separations. From the mock catalogues, we find 
that the fraction of random projected pairs with separations 
< 15 h-^kpc averaged over the 26 fields is 0.0002 ± 0.00009. 
As a result, the SSCpairs are not random projections. 

To estimate the interloper contribution to the pair sep- 
arations, we minimised a of ^ two component dis- 
tribution that includes the random catalogue results and 
the angular clustering of LBGs. For LBG angular cluster- 
ing, we used the conventional form of the angular correla- 
tion function, uj{0) — A6~^ . We first test the fit by fixing 
the slope of the angular correlation function to 5 = 0.57 as 
found by lAdelberger et al. (2005) and varying the fraction 
of interlopers, //. We find the best fit when using an in- 
terloper fraction of // = 0.196 and // = 0.195 for the C05 
and S03 samples, respectively. The contribution by the LBG 
angular correlation function is shown in Figures [6] and [71 as 
the thick dashed (light blue) lines with the combined two- 
component LBG/random contribution (corrected angular 
correlation function) shown by the thick solid (blue) lines. 
To help illustrate the effect the slope of the angular corre- 
lation function has on the results, we also show the LBG 
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Figure 7. Same as Figure[6]but in log/log space to better discern 
the behaviour on small (< 20 h~^kpc) scales. The photometric 
close pairs (solid histogram) fail at separations smaller than ~ 7 
h~^'kpc because of the resolution of the images (also see text). 
The number of close pairs with the addition of the SSCpairs is 
represented by the dashed histogram. The top panel shows the 
close pairs and SSCpairs presented in thi s work. The c e ntral panel 
shows the close pairs for the survey of ISteidel etHI (|2003l i and 
includes predicted SSCpairs using an expected distribution based 
on the observed distribution found in our survey. The bottom 
panel shows the combined datasets. These plots illustrate that the 
predicted SSCpairs fall within the expectations of all functional 
fits at small radii and that the probability that the SSCpairs are 
the result of random pairings is extremely low. 
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contribution when using a slope of ^ = 0.8 (dashed green 
Une) . 

We re-fit the data letting both parameters vary freely 
and found best fits in 6 and // within ^ 3% of the 
above values. Because this type of test has a small de- 
pendence on choice of bin size, we use the fixed S = 0.57 
(jAdelberger et al]|2005l l result in our assessment of SSCpairs 
expectations. The consistency in pair distributions and pa- 
rameter fits for both surveys illustrates the utility of using 
pair separations to describe the behaviour of LBGs on small 
scales. 

The number of LEG close pairs predicted by the cor- 
rected angular correlation function provides a good match 
to the observed pairs on all scales down to the resolution 
of the images and becomes indistinguishable from the un- 
corrected LEG angular correlation function expectations at 
separations smaller than ^ 30 h~^kpc. This agreement reit- 
erates the negligible contribution of random projected pairs 
on small scales. 

The widths of our spectroscopic slitlets are < 8 /i~^kpc. 
Inspecting Figure [71 at separations of < 15 h~^kpc shows 
that the expectation from every functional fit to the data 
predicts that we should find close pairs serendipitously in 
our slitlets. From the fit of the corrected (and uncorrected) 
LEG angular correlation function, we find that both the ex- 
pected number of SSCpairs and the expected distribution of 
SSCpair separations agree well with that of the observations 
(shown as the dashed histogram in Figure [7]) . 



5 SERENDIPITOUS PAIR EXPECTATIONS 
FROM SIMULATION 

We test the identification of the SSCpairs against the predic- 
tions of LCDM by analysing a hi^h-resol ution hybrid numer - 
ical/analytical cosmological simulation teer rier et al.|[2006l V 
The full description of the simulation a nd the close pa i r ana l- 
ysis is presented in a companion paper jEerrier et al.1 (|2009l l. 
Eelow we describe a few of the relevant specifics and results. 

For large-scale structure, th e simulation uses an adap- 
tive refinement tree N-body code (Kravtsov et al. 1997) that 
follows the evolution of 512^ particles in a comoving box 
120 Mpc on a side. Numerical overmerging on small 
scales is overcome b y incorporating the analytic method of 
Zentner et al. to trace the substructure. We sample 

the simulation at the z — 3 time-step and use dark mat- 
ter halos in the simulation that match the spatial corre- 
lation function of the form ^(r) = (r/ro)~^ measured for 
the uiR ^25.5 LEG samples. We adopt the spatial correla- 
tion length ro = 3.3 ± 0.6 ([Coo ke et al. 2006) using a fixed 
slope 7 = 1.^ f or the C05 survey and ro = 4.0 ± 0.6 and 
7 = 1.57 ± 0.14 flAdelberger et al.ll2005h for the S03 survey. 
We find that halos with Vdrc ^ 140 km s~^ in the simula- 
tion yield a correlation length of ro = 3.9 it 1.5 and slope 
of 7 = 1.6 ± 0.3. These halos correspond to a number den- 
sity of riLBG — 0.016 Mpc~^, comoving, and halo mass 

4 In lCooke et HI (l2006l ) we use several methods to measure and 
test the correlation function parameters because of the inherent 
uncertainties in analysing ~ 200 LBGs. We choose the value we 
find most reliable for the 140 LBGs used here that is obtained 
from a maximum likelihood analysis holding 7 fixed. 



Muaio > lO^^-^Mo {{MHalo) - lO^^'^M©) which are in very 
good agreement with the corrected volumetric density from 
the photometric sample and the mass inferre d from clus- 
tering statistics lO^^-^Mrr.. Adelberger et aLll2003l . liooHl : 
ICooke et al.ll2006h . 

The simulation is then analysed in a manner identical 
to the observations using (1) the constraints from the see- 
ing FWHM and centroid accuracy from the images, (2) the 
MOS spectroscopic slitlet dimensions of our survey, (3) the 
spectral resolution, and (4) the uncertainties in determining 
the systemic LEG spectroscopic redshifts. We compute the 
number of pairs that fall in mock slitlets with random orien- 
tations that are resolvable in the images and that have ve- 
locity offsets that would be detectable in LRIS spectroscopy. 
We find an expected fraction of 0.052 ± 0.002, or an aver- 
age of 7.1 LEGs (3 — 4 pairs) that is to be compared to the 
0.069 ±0.023 fraction found from the analysis of this survey. 
This close agreement helps reinforce the interpretation that 
the SSCpairs are indeed close z ^ 3 LEG pairs as a similar 
number are expected from LCDM to be found serendipi- 
tously in our survey. We find that the expected fraction of 
SSCpairs from the simulation is rather insensitive (changes 
by a few percent) to a variation of the LEG correlation pa- 
rameters within the range of their uncertainties. 



6 SERENDIPITOUS LY-a EMITTERS 

Lya emitters (LAEs) are defined here as objects display- 
ing detectable Lya emission but have continua fainter than 
R — 25.5 and therefore are not candidates for inclusion in 
the LEG spectroscopic sample. We discovered two serendip- 
itous z - 3 LAEs in the MOS slitlets of the fuh C05 LEG 
sample. This is approximately the number we would expect 
to fall serendipitously into our slitlets, when we consider (1) 
the depth of our spectroscopy, (2) the ~ 1 arcmin^ solid 
angle subtended by our slitlets, and (3) the z ^ 3 LAE 
densities reported by ICawiser et al.l (|2007[ ) and I Grove et al.l 
(|2009l ) when assumed to be constant over the redshift path 
probed. The detected objects are LAE-228 at 2; = 3.007 in 
field PC1643+4631A and LAE-475 at z = 3.209 in the field 
ERI1013+0035. Figure [8] presents their 1-D and 2-D spec- 
tra. Eoth redshifts are measured from their Lya emission 
peaks and exhibit Lya EW of 48 A and 140 A, respectively. 

Interestingly, LAE-228 was found as a spectroscopic 
close pair to the spectroscopically confirmed LEG-228 with 
z = 3.008. The projected separation between the two sys- 
tems is 22.7 kpc and they exhibit a AvLya = 450 km s~^. In 
addition, LAE-475 was found near LEG-475 with z = 3.218 
having a projected separation of 40.1 kpc and a Lya veloc- 
ity offset of 650 km s~^. Although both systems could be 
considered close pairs, neither is included in the analysis in 
this paper because the LAEs do not meet the uir ^ 25.5 
criterion and becau se LEG-475 is con sidered a "probable" 
LEG (category 4 in lCooke et"alll2006l ) from the lower S/N 
and quality of the spectral-line fit. 

LAEs are believed to be systems with lower mass and 
lower star formation rates than that of the LEG popula- 
tion (- 1O^°-^~^^M0 and ^ lOM© yr"^; e.g.. ICawiser et al.1 
I2OO7I ) that are undergoing a recent star forming episode. The 
fact that the two LAEs detected are close pairs to one, and 
likely both, LEGs suggests that a fraction of the Lya emis- 
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Figure 8. Spectra of the LAEs and close LBGs detected in the 
survey. Top split panel: 1-D spectra of LBG-228 and LAE-228 
showing the Lyo; emission and continua of both objects. Cen- 
tre split panel: 1-D spectra of LBG-475 and LAE-475. LAE- 
475 show strong Lyo; emission, however, LBG-475 is considered a 
probable LBG because of the low S/N and inconclusive nature of 
the ISM line fit. Bottom two panels: The 2-D spectra for pairs 
228 and 475. The position of the Lyo; features and continua for 
pair 228 and 475 are indicated and have projected separations of 
22.7 and 40.1 h~^kpc, respectively. Although they are within the 
close pair separation criteria, they are not included in the analysis 
because the LAEs do not meet the m^^ < 25.5 criterion. 



sion, and therefore the detect ability of LAEs, may result 
from interactions. 



7 Lya EMISSION AND CLOSE PAIRS 

In our analysis, we see dominant Lya in emission from all 
galaxies of the SSCpairs and a prevalence for dominant Lya 
emission for LBGs having close photometric companions 
(SSCpair candidates; see Figures [1] & [3]) . We consider the 
frequency of the appearance of Lya in emission for the SSC- 
pairs and the SSCpair candidates. The fraction of eLBGs in 
our spectroscopic sample of 140 LBGs is 0.55. This fraction 
includes LBGs exhibiting purely emission with undetectable 
absorption, to those with a complex profile that show Lya 
absorption but have strong enough emission to result in net 
Lya emission EW. By our definition, the larger spectro- 
scopic sample of SOS finds a comparable eLBG fraction of 
0.51. We remark that the low S/N, low-resolution LBG con- 
tinua and complex Lya features near zero net EW, where 
the peak of the EW distribu tion is located , may produce 
minor differences between the lShaplev et al 

] (Hqos) analysis 

of the SOS data and ours. These diflFerences may contribute 
toward the small discrepancy in the two eLBG/aLBG ra- 
tios, however, do not affect the results below for LBGs with 
< 15 /i~^kpc separations as none of these exhibit such Lya 
profiles. In the probability calculations below, we use the 
eLBG fractio n from our survey exc lusively. The eLBG frac- 
tion found by Shap lev et al. analysis of the SOS data 
would result in lower probabilities. 

The binomial probability of observing all 10 SSCpair 
galaxies with dominant Lya emission and zero galaxies with 
dominant Lya in absorption by chance is p = 0.0025. The six 
spectroscopic LBG components of the SSCpair candidates 
all show dominant Lya in emission but were not included 
in the SSCpairs analysis because the bulk of the secondary 
component did not enter the slitlet (Figure [1] and Figure [S). 
Combined, the binomial probability that all 16 galaxies are 
eLBGs and zero are aLBGs is p = 7.0 X 10"^ If the emis- 
sion of these systems was a result of star forming clumps or 
single LBGs, we would expect to observe an approximately 
equal fraction (or ^ 45% for this survey) of regions that are 
dominated by Lya absorption and we observe none. As a 
reminder, systems exhibiting Lya absorption (aLBGs) were 
searched for thoroughly in the data f ^S.2|) . 

When examining the Lya features of spectroscopic 
LBGs having photometric companions that meet our colour 
criteria and have projected separations between 15 — 50 
/i~^kpc, with the caveat that some of these may not be 
physical pairs, we find 11 eLBGs and 7 aLBGs {p = 0.17). In 
addition, if we include all spectroscopically confirmed LBGs 
with photometric close pairs that have relaxed colour crite- 
ria (see C05) and that have projected separations between 
15 — 50 h~^kpc (those with < 15 h~^kpc are included in 
the SSCpair candidates), we find a total of 15 eLBGs and 8 
aLBGs (p = 0.11). These two measurements suggest an over- 
abundance of eLBGs and, combined with the abundance of 
Lya emission in the SSCpairs and SSCpair candidates, indi- 
cate an overabundance that diminishes with increasing ra- 
dius. Figure [9] plots the number of eLBGs and aLBGs with 
photometric companions with separation. Finally, LBG-228 
with close pair LAE-228 (22.7 h~^kpc) exhibits Lya in emis- 
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sion. The tentative LBG-475 with close pair LAE-475 (40.1 
h~^kpc) shows weak evidence of complex Lya absorption 
and emission in the low S/N spectrum, but a proper EW 
analysis is not possible at this time. 

We posit that the overabundance of eLBGs at small 
separations (< 15 h~^kpc) is caused, in part, by induced 
star formation and the dispersal of gas and dust rela t ed in - 
teractions. Object D3 in the sample of Shapley et al. (2003|) 
provides further support in that both galaxies in the iden- 
tified close pair (~ 15 h~^kpc separation) exhibit dominant 
Lya emission. The extended Lya emission of D3, Lya blob 
nature, and potential fraction of escaping Lyman continuum 
photons reinforces this scenario. Finally, the trends of Lya 
emission with morphological parameters are consistent with 
a picture where a non- negligible fraction of LBGs produce 
double-peak Lya emission when merging their ground-based 
1-D spectra as a result of triggered or more revealed Lya 
emission. 

The statistics presented here can be quickly improved 
from the analysis of other similar existing and future 
datasets. Nevertheless, the Lya behaviour of the SSCpairs 
appears real. Such behaviour would suggest that although 
star formation is ubiquitous at high redshift, caused in part 
by the high gas fraction in most systems, and generates Lya 
emission, close pairs have additional Lya emission and UV 
flux as a result of an interaction. The enhanced Lya emis- 
sion from interactions modifies the ratio of eLBGs and in- 
creases the fraction of detectable eLBGs for a given criteria 
whereas without an interaction they would remain either 
too faint or too red. Such behaviour would bias any colour- 
selection technique. Similarly, interactions may make LAEs 
detectable whereas without an interaction they would re- 
main too faint to be detected from their continuum alone. 

We find ~ 7% of the galaxies in our spectroscopic sam- 
ple identified as interacting pairs. If the higher fraction of 
Lya double-pe ak systems in t he hi gher S/N, higher resolu- 
tion sample of IShaplev et"aD (|2006h prove to be interacting 
pairs, because both samples are serendipitous detections, the 
true fraction of interacting LBGs at z 3 may be > 30%. If 
borne out, the double Lya emission/double morphology ob- 
servable in ground-based data may provide a simple means 
to identify a significant fraction of interacting galaxies and 
help to constrain the merger rate at high redshift. The rela- 
tionship between Lya emission and pair separation will be 
fully explored in a future paper. The evidence is presented 
here to establish a link between the presence of Lya in emis- 
sion and the identification of interacting pairs. 

In conclusion, the above picture and the observations of 
merged 1-D double Lya emission profiles in ground-based 
data as markers of this process are supported by: (1) the 
detection of Lya in emission from every component of the 
SSCpairs and the spectroscopic components of all six SSC- 
pair candidates, (2) an indication of an overabundance of 
eLBGs in spectro-photometric ^ 50 h~^kpc close pairs, (3) 
the trends of Lya emission EW with nucleat i on, co lour, and 
separation (size/multiplicity) of iLaw et al.l (|2007f ), (4) the 
plausibility from Lya and ISM feature analysis that the dou- 
ble peak systems in the work of Shapley et al. (2006) are 
close and/or interacting pairs, (5) the fact that both galax- 
ies in the identified close pair D3 exhibit Lya in emission 
and the identification of this system as an extended Lya 
blob, (6) the potential Lyman continuum escaping photons 
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Figure 9. Histograms of the Lyo; behaviour of spectroscopic 
LBGs with close photometric companions. Pair separations where 
the spectroscopic component exhibits dominant Lya absorption 
indicated by the red (grey) back-hatched histogram (upper panel) 
and dominant Lyo; emission are indicated by the blue (black) 
forward-hatched histogram (lower panel). The SSCpairs are not 
included. All SSCpair components have separations of < 15 
/i~-'^kpc and show Lyo; emission. Pairs shown with < 15 /i~-'^kpc 
separations include the SSCpair candidates. 



from D3 and double-peak Lya emission system C49, and (7) 
the detection of one, and potentially two, LAEs within ^ 50 
/i~^kpc of spectroscopic LBGs in this survey. 



8 SUMMARY 

A systematic search of deep Keck imaging and spectroscopy 
of 140 z ~ 3 LBGs in nine separate fields has uncovered five 
serendipitous spectroscopic close pairs (SSCpairs) with sep- 
arations of < 15 /i~^kpc. Below, we summarise our findings 
and the implications of the imaging and spectroscopic data 
from this survey, the photometric data from the survey of 
S03, and the results from an analysis of a carefully matched 
high-resolution cosmological hybrid simulation. 

(1) We find evidence in the imaging and spectroscopic ob- 
servations that the SSCpairs are two distinct and likely in- 
teracting systems. 

(i) The SSCpairs exhibit two fiux peaks with ^3 — 13 
/i~^kpc projected separations in the broadband images. The 
VRI filters probe the continua of O and B stars at z 3 and 
the SSCpair components exhibit luminosities that are typi- 
cal of z ^ 3 LBGs. In addition, morphology of the SSCpairs 
is consistent with the expectations of low to intermediate 
redshift UVLGs projected to z ^ 3. UVLGs are in every 
way analogues to LBGs having similar mass, gas content. 
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UV colours, metallicity, half-light radii, and star formation 
rates. Mergers/interactions are determined to be the main 
star- forming mechanism for UVLGs in the sample of O08. 

(ii) Each component of the SSCpairs exhibits detectable 
Lya emission and a corresponding 2-D continuum. The con- 
tinua of the SSCpairs are closely spaced or overlapping. 
When merging the ID spectra, each SSCpair shows a double- 
peak Lya emission profile which is therefore different in na- 
ture than those modeled for a static or expanding shell of gas 
driven by galactic-scale outflows. The Lya emission peaks 
and continua in the 2D spectra are spatially offset and ex- 
hibit a one-to-one spatial correspondence with continuum 
flux peaks in the images. 

(iii) The velocity offsets between the two Lya emission 
peaks of the SSCpairs (before and after angular separation 
corrections) are not representative of a random sampling 
of Lya offsets pulled from the LBG population if the two 
features are generated by two sources within single systems. 
Instead, the Lya offsets are in agreement with single Lya 
emission from two separate systems. 

(iv) We see evidence for two sets of ISM absorption lines 
(at a lower confidence level) in each merged SSCpairs spec- 
trum having relative velocities inconsistent with the expec- 
tations of the general LBG population but consistent with 
two separate systems. In addition, the velocity differences 
between the Lya and ISM features assigned to each compo- 
nent of the SSCpairs, under the assumption that they are 
distinct systems, is consistent with typical LBG values. 

(v) The high-S /N, higher-res o lution 2; ~ 3 LBG spectro- 
scopic sample of IShaplev et al.l (|2006l l report that four of 
the nine eLBGs in their sample of 14 systems display double- 
peak Lya emission. Analysis of the reported values finds the 
data for the three double-peak systems and one identified 
double system to be consistent with the SSCpairs observa- 
tions. Both samples show Lya and potential ISM offsets and 
Lya -ISM velocity differences that are less consistent with 
a single system and more consistent with typical values for 
two interacting z ^ 3 galaxies. 

(vi) The morphology and Lya features of the SSC- 
pairs are consis tent with the reported relationships of 
I Law et aD (|2007l l between space-based morphological pa- 
rameters (Gini, multiplicity, size, and colour) and Lya emis- 
sion EW of z ~ 3 LBGs. An interpretation of the these rela- 
tionships points to a picture where single and double-nuclei 
LBGs are Lya emitting LBGs, with a significant fraction 
being interacting systems. 

(2) We examine the distribution of all close pairs in the data 
to determine the expectation of serendipitous pairs in our 
spectroscopic slit lets. We find that an angular correlation 
function of the form uj{0) = A6~^ with slope 6 = 0.57 pro- 
vides an excellent fit to close pair counts for both surveys 
tested. The fit incorporates contamination by random in- 
terlopers determined from mock catalogues matched to the 
volume and behaviour of the data. The slope of the fit to 
the close pair counts is identical to, or in very close agree- 
ment with, the z ^ 3 LBG values in the literature derived 
from angular and spatial clustering statistics. Extrapolation 
of the fit to below the resolution limit of the photometric 
data (< 10 h~^kpc) predicts the existence of SSCpairs with 
a similar number and separation distribution to that ob- 
served. 



(3) We analyse a high-resolution cosmological hybrid nu- 
merical and analytical simulation to assess the expectation 
of the SSCpairs from the predictions of LCDM. We care- 
fully matched the simulation to the number density, spatial 
correlation function, spectral and imaging resolutions of the 
data, and the dimensions of the slitlets. The analysis finds 
a prediction of ^ 4 pairs to fall serendipitously in the multi- 
object sUtlets of our survey, in good agreement with the 
observations. 

(4) We present two LAEs discovered serendipitously in the 
MOS spectroscopic slitlets which are expected from LAE 
densities at ^ ~ 3 and the solid angle subtended by our 
slitlets. However, both LAEs are ^ 50 h~^kpc close pro- 
jected pairs to LBGs. We confirm one of the LAEs to lay at 
the redshift of a nearby LBG and is a bona-fide close pair 
(22.7 /i~^kpc), whereas the second is only a probable close 
pair because of the low S/N of the LBG spectrum. In ad- 
dition, the LBG in the bona-fide close pair exhibits Lya in 
emission. 

(5) We uncover a relationship between Lya emission and 
pair separation. We find an overabundance of LBGs with 
separations of < 50 h~^kpc exhibiting Lya in emission as 
compared to the full sample. Moreover, without exception, 
all galaxies in the SSCpairs and all SSCpair candidates 

15 h~^kpc ) have dominant Lya emission. 

The presence of Lya in all SSCpairs and SSCpairs can- 
didates supports the interpretation that the SSCpairs are in- 
teracting LBGs under the assumption that interaction may 
disperse gas and dust and trigger observable star formation. 
Because of the close proximity of the galaxies in the data 
to each other, the double Lya emission/double morphology 
of the SSCpairs provides a useful and consistent diagnostic 
to help detect interacting close pairs from 2 < 2; < 7 in 
similarly obtained existing and future optical ground-based 
data. 
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